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Abstract 

We apply linear-response analysis of the Gross-Pitaevskii equation to ob- 
tain the excitation frequencies of a Bose-Einstein condensate confined in a 
time-averaged orbiting potential trap. Our calculated values are in excellent 

agreement with those observed in a recent experiment. 
PACS Numbers: 3.75.F1, 67.40.Db, 67.90.+Z 
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The recent attainment of quantum degeneracy conditions in magnetically trapped alkali 
vapors [|1],0|§ has opened the road to understanding of the many-body physics of Bose- 
Einstein condensates (BECs) in unprecedented detail. For dilute gases, it is believed that 
the essential physics of the BEC ground state is captured in the Gross-Pitaevski mean-field 



agreed reasonably well with the few experimental determinations of condensate shapes, 
sizes, and lifetimes that have been made to date, but it cannot be said that the theory has 
been subject to stringent tests. 

In this paper we report theoretical results for the excitation spectrum of a BEC of 
trapped 87 Rb, obtained by computing the response to small mechanical disturbances of a 
BEC described by the GP equation ||. These results are compared with those of a recent 
experiment [0], which has observed the free oscillations of a BEC that is briefly shaken 
at frequencies near resonance. We believe that this comparison provides the most critical 
quantitative test of mean-field theory made to date. The agreement between experimental 
and theoretical results is excellent. This suggests that the GP equation and its variants can 
adequately describe excited-state as well as ground-state properties, and so should provide 
a practical framework for explorations of emerging issues of condensate dynamics, finite- 
temperature phenomena, and atom laser || design. 

To describe a magnetically trapped atomic gas, we adopt the standard GP equation, 
which is applicable @.|1(J when the condensate fraction of a gaseous system is close to 
unity. Each atom in the condensate occupies the same orbital ip g {r), which is determined 
by solution of the nonlinear Schrodinger equation, 



where H = — |j^V 2 + Vtra P (r) is the Hamiltonian for an isolated atom in the trap, iVo is 
the number of atoms in the condensate, Uq represents the interaction between condensate 
atoms, and the eigenvalue \x is the chemical potential. 

In most current trap designs, Vt ra p can be described by the anisotropic harmonic oscillator 



formalism. Calculations done with the Gross-Pitaevski (GP) equation have indeed 
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potential, Vt rap (r) = ~m (lo^x 2 + uj 2 y 2 + uj 2 z 2 ^j, where m is the atomic mass and a>i = 2iu>i 
is the angular frequency of oscillation along the axis i. The time-averaged orbiting potential 
(TOP) trap [[[jn] treated here is cylindrically symmetric; its potential is given by uo x = uo y = 
lu_l and uj z = \^8tu±. The results presented below correspond to v z = lu z /2tt = 210 Hz. The 
parameter Uq expresses the interaction between two atoms as Uq = 47ih 2 a/m, where a is the 
scattering length, which characterizes the zero-energy behavior of the s-wave phase shift in 
collisions between two atoms. The scattering length a is the only piece of atomic collision 
data used as input to our calculations. The present results are given in terms of the most 
recent ]15|] experimental value, a = llOao, where a is the Bohr radius; our calculations 



were actually carried out with a previously published [0 value of a = 100a , but since 
Eq. m obeys a scaling law involving iVo, u±, and a (see below), we can rescale our results 
to compare quantitatively with experiment. It should be noted that for the alkali atoms in 
current BEC studies, the experimental determination of a requires extensive spectroscopic 



analysis, and present values are accompanied by substantial uncertainties |T2 . 

Equation [I] has previously been solved by several independent methods [|],[5] to describe 
the BEC ground state. Here we investigate the response of the ground state to an oscillatory 
perturbation at angular frequency u p . The associated time-dependent GP equation takes 
the form: 

ih™ = [Fo + f/o|^(r,t)| 2 + / + (r)e-^ + /_(r) e ^]M/(r,t), (2) 

where /±(r) are the spatially-dependent amplitudes of the perturbation. We solve this 
equation in the linear-response limit. The details of this approach are described elsewhere 
|13fl, and we simply state the central results here. By using the Ansatz 



*(r,t) = e 



-i/it/H 



(3) 



N^ g (r) + u{r)e' iUpt + v{r)e luJpt 
we obtain the linear-response equations, 

\C - huj p ] u(r) + N U [^(r)f v(v) = -/+(r)^(r), (4) 
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N U r g (r) u(r) + [£ - fto;,)] v(r) = -/_(r)^(r), 



(5) 



where £ = H - ^ + 2U N Q |^(r)| 2 . 

This pair of equations can be solved by expansion in terms of the GP normal-mode 
equations, 



[C - hco x ) u x (r) + NqUq [^ 9 (r)]% A (r) = 0, 



(6) 



N U 



C(r) u A (r) + [£ + hw x ] v x {r) = 0, 



(7) 



where uj\ is an eigenvalue and u\(r), V\(r) are corresponding eigenfunctions. For a homo- 
geneous BEC (i.e. V^ rap (r) = 0), uj\ would be a quasiparticle excitation frequency with 
a continuous spectrum, corresponding to travelling sound waves. In a trap, however, the 
elementary excitations of the BEC remained confined, so u\(r) and v\(r) must be square- 
integrable, and the spectrum of u\ is thus discrete. The solution of eq. |3| can then be 
obtained in the ususal way by a superposition of eigenfunctions: 



u[r) 
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t ux(x) ^ 



(8) 



where 



g x = / d 3 r [ul(v)f + (v) +<(r)/_(r)]V s (r). 



(9) 

The response diverges when uo p = u\, i.e. when the driving frequency is equal to a natural 
excitation frequency. In practice, of course, this divergence will be eliminated by damping 
effects, a subject of ongoing work that will not be discussed here. To the extent that mecha- 
nisms of damping and their associated frequency shifts can be neglected, it is apparent that 
free oscillation of the condensate after a transient disturbance will occur at the frequencies 

As we have discussed elsewhere [T3], there is a straightforward connection between the 
resonant oscillation frequencies uj\ and the the quasi-particle mode frequencies that are 



encountered in the standard Bogoliubov approximation []T0| . Stated simply, Eqs. (|) and (^) 



are identical to the equations that define the quasi-particle modes and frequencies within the 
Bogoliubov approximation. Thus an experiment that measures the free oscillatory response 
of a shaken BEC provides a direct observation of the quasi-particle spectrum. In particular, 
by shaking the BEC at a frequency near one of the resonances u)\, one can produce a response 
that is dominated by the u>\. This is the approach that has been taken by the first such 
experiment reported 0. 

We have solved numerically the system of equations consisting of Eqs. ([!]), (H), and 
(0) under the conditions of that experiment. The solution was accomplished in two steps. 
First, Eq. (HD was solved by expanding the solution ijjg{r) in a basis set consisting of a finite 
number of trap eigenfunctions. The details of the numerical method have been recounted 
elsewhere ||. Equations (||), and (0) were then solved by expanding w A (r) and v\(r) in 
the same basis set. These expansions convert Eqs. (P), and (0) into a generalized matrix 
eigenvalue problem that can be solved by standard numerical techniques The error 

in the solutions was assessed by increasing the basis-set size until the mode frequencies 
converged to at least three figures. 

Figure |I| shows our results for the lowest three excitation frequencies (in units of the 
trap frequency v±), as a function of Nq. A simple scaling law faciltiates comparison of 
calculation and experiment. A solution of Eqs. ([!]), (H), and ([?]), plus normalization for 
experimental values of the parameters {N^ e \ a e , m e } } will also satisfy the equations for 
the parameter set {Nq , a t , , m t }, if the quantity 7 = N a (mz/jj 1 / 2 is constant. Thus, 
excitation measurements performed on a BEC with iV (e) atoms in a trap of frequency 
can be related to the spectrum displayed in Fig. [1] by taking 

(e)\V2 

ATM 



= do) 

where vf = 210/ « 74.25 Hz. 

We have used this scaling law to compare our results with those of the recent experiment 
where excitations of an atomic BEC were observed for the first time. Table | presents 



the comparison. No attempt has been made here to account for experimental uncertainties. 
The ranges presented for the experimental ratio fz^de/^.! ) involve a 1-3 % reduction of 



the raw data to extrapolate to zero-amplitude driving force |14[]. It is fair to say that the 
agreement is excellent, as the difference between theory and experiment ranges from 2-6 %. 

To better understand the nature of these excitations we compare them with the results 
of Stringari 0] who obtained analytic solutions to the linearized GP equation in the hy- 



drodynamic limit (iVo — > oo). The curve labeled "m = 1" is a doubly-degenerate dipole 
excitation that coincides exactly with the first excited state of the bare trap. This is because 
the lowest dipole mode of an ensemble of identical interacting atoms in an external harmonic 
potential corresponds to a rigid motion of the center-of-mass, independent of the nature of 
interatomic forces [17] . 



The curve labeled "m = 2" corresponds to doubly degenerate excitations that correlate to 
the lowest m = |2| energies of the two-dimensional harmonic oscillator in the noninteracting 
limit. The quadrupolar nature of these excitations is exhibited in fig. |2], which contains a 
plot of u\ over a region of the xy plane [i.e. u\(x, y,0)} for N = 2000 atoms. This plot 
clearly shows the four-peak structure characteristic of quadrupole excitations. The large- N 
limit of this mode, is m ode/ v± — > y/2, is also shown in figure [I] The fact that the middle curves 
also seem to be approaching this limit as iVo — > oo also lends support to their interpretation 
as quadrupole excitations. 

The curve labeled "m=0" corresponds to a "large-iVo" excitation that is a breathing 

mode in the xy plane. The asymptotic limit is shown in fig. p] and is given by ]17| 

/ 3 i , \i/2 

iwK - (2 + -A 2 - -V9A 4 - 16A 2 + 16} (11) 

where A = v z fv± = \/8 for the TOP trap. As discussed in ref. 0, the symmetries of the 
normal modes can be tested by experimental selection rules, and the classifications of the 
observed modes are found to agreement with those given here. 

In conclusion, we have presented excitation spectra that agree well with the data of a 
recent experiment. We have shown that these data constitute a direct measurement of the 
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T = Bogoliubov spectrum of an atomic BEC. We have also used the large- N limit and 
mode shapes to describe the nature of these excitations. The experimental confirmation of 
these data will have significant implications for understanding of the many-body physics 
of these dilute, weakly-interacting bosonic systems, and for practical use in future BEC 
engineering. 
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Figure 1, "Collective excitations of atomic Bose-Eiiislein condensates", 
M. Edwards, R.J. Dodd, C.W. Clark, P.A. Ruprechl, and K. Burnett 



FIG. 1. The lowest three calculated excitation frequencies, in units of the perpendicular trap 
frequency Vj^ 7 of the JILA TOP-trap condensate, it vs. the number of condensate atoms, Nq\ 




Figure 2, "Collective excitations of atomic Bose-Einstein condensates," 
M. Edwards, R.J. Dodd, C.W. Clark, P.A. Ruprecht, and K. Burnett 

FIG. 2. The UA~component of the m = |2| excitation (middle curve of fig. ||), displaying its 
quadrupolar shape. 
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TABLES 

TABLE I. This table presents the comparison with excitation data obtained in the experiment 
of ref . [4] . 



mode 


(e) 


(e) 

mode 








(t) 
v j 

nioac 
* 


% diff. 


m = 


43.2 


82 


3420 


1.84—1.88 


2870 


1.89 


2.7 


m = 


132 


244.7 


2400 


1.79—1.83 


3520 


1.88 


5.0 


m = 2 


43.2 


62.5 


2800 


1.41—1.44 


2350 


1.49 


5.7 


m = 2 


132 


188.5 


2200 


1.39—1.42 


3230 


1.47 


5.8 
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